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I. INTRODUCTION 



Since observations show that space is homogeneous at large scales to a very good approx- 
imation, homogeneous models have been studied for many decades in classical and quantum 
cosmology They provide insights into the behavior of our universe while avoiding all 
the complicated field theoretic details of full gravity by their restriction to finitely many 
degrees of freedom. Despite the huge reduction by infinitely many degrees of freedom, a 
large variety of different models is left which allow a detailed investigation of various non- 
trivial issues such as: Dirac observables, the embedding of symmetric models (e.g., isotropic) 
into less symmetric ones, dynamics, the approach to classical singularities and cosmological 
phenomenology. 

Loop quantum cosmological methods for those models have recently been developed 
and shown to simplify considerably after a diagonalization of the connection and triad 
degrees of freedom jfj. As an example the dynamics of the Bianchi I model has been studied 
and shown to be singularity-free. In this case the Hamiltonian constraint resembles that of 
the full theory |j| (albeit it is much simpler), and thus the quantization can be regarded 
as a reliable test. It turned out that the extension from the isotropic case j^, 0] to this 
anisotropic case, requires certain features which are in fact present in quantum geometry. 
In particular, the mechanism of a singularity-free evolution in anisotropic models relies on 
the fact that the evolution extends through the classical singularity which can be identified 
with a submanifold in minisuperspace. This requires the classical singularity to lie in the 
interior of minisuperspace, rather than at the boundary (including infinity). Here it is 
essential that quantum geometry is based on densitized triad variables where in fact the 
Kasner singularity is in the interior (all densitized triad components being zero), while it 
would be on the infinite boundary in co-triad variables. The isotropic case is not sensitive 
to this issue since both densitized and co-triad variables have the classical singularity in 
the interior of minisuperspace (here, it is essential to use any (co-)triad rather than metric 
variables); thus the existence of the extension of the methods to a homogeneous model is 
non-trivial. An extension to the full theory of the general class of inhomogeneous models 
seems much more complicated at this point. However, it has been conjectured that at 
the classical level, the most general homogeneous behavior is that of the Bianchi-IX model 
and furthermore that it describes the approach to the singularities even of inhomogeneous 
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models. It is therefore important to study the more general Bianchi class A models. 

Homogeneous models other than Bianchi I, however, present an additional complication 
since they have non-zero intrinsic curvature and, as a consequence, their spin connection 
cannot be zero. This is in contrast even to the full theory where the spin connection can 
be made arbitrarily small locally by choosing appropriate coordinates. In a homogeneous 
model the freedom of choosing coordinates is restricted to those which preserve (manifest) 
homogeneity, and this implies that the spin connection is a covariant object and that it has 
to have a certain size in a given model. This has to be taken into account properly when one 
quantizes the Hamiltonian constraint. In particular, in a classical regime only the extrinsic 
curvature can be assumed to be small, but not necessarily the intrinsic curvature which 
determines the spin connection. This is essential to understand the semiclassical behavior. 

Since a special treatment is required which is not necessary in the full theory, one has to 
be more careful when interpreting the results. It introduces more quantization ambiguities 
which have to be shown not to influence main results. On the other hand, more possibilities 
for phenomenology emerge which, at least qualitatively, can often be seen to be insensitive 
to ambiguities. 

Having a more distant relation to the full theory, it is helpful to have a set of admissibility 
criteria. One such criterion is of course to accommodate a semiclassical approximation to 
the quantum dynamics. Since a semiclassical description is based on continuum geometry 
while the quantum dynamics is in terms of discrete quantum geometry, this criterion is 
formulated here in terms of a continuum approximation. A further criterion is provided by 
the requirement of a locally stable evolution 



evolution equation derived along the lines of 



1|. Such a requirement arises because the 
is usually a high order difference equation 



which at large volume is well approximated by the second order Wheeler-DeWitt differential 
equation. Hence, one can always construct solutions by choosing initial data, at large volume, 
to be that provided by solutions of the Wheeler-DeWitt equation. Perturbation of such an 
initial data will generate solutions which will also include surplus solutions of the high order 
difference equation. Generically these will have Planck scale oscillations. There is then 
the possibility that such extra solutions can become dominant under the evolution. In 
fact these are expected to become dominant as one gets to smaller volume since there are 
huge differences between the continuum and the discrete formulations in the Planck regime. 
However as one evolves to larger volumes, the perturbations must not grow too rapidly 
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since this would imply domination by solutions with Planck scale oscillation even in the 
classical regime. The requirement we are looking for, local stability, prohibits this behavior 
by demanding that the local behavior of solutions to the difference equation around a large 
value of the evolution parameter is not exponentially increasing. In this paper we present 
a quantization of the Hamiltonian constraint for Bianchi models with non-zero intrinsic 
curvature which fulfills this condition and has the correct semiclassical limit. The first 
quantization given in ^ was not admissible in this sense; thus, the quantization given there 
is valid only for the Bianchi I model (for which we have the same quantization here). 

It turns out that the local stability condition is selective: it requires all roots of a high 
order polynomial to have unit norm which is not easy to achieve randomly. The selectivity is 
increased by the fact that the same strategy of quantizing a Hamiltonian constraint should 
work in all homogeneous models, using different procedures in different models would imply 
that not all of them can be related to the full theory. 

The stability condition can be side-stepped by quantizing connection components in the 
constraint by hand such that only a second order difference equation results. While this 
would eliminate local stability as a selection criterion, it would also imply that the quanti- 
zation is even more distant to that in the full theory where such a quantization cannot be 
possible. Further clues as to the necessary order of the evolution equation can come from 
studying Dirac observables, which will be pursued elsewhere. 

In section |Hj we briefly describe the diagonal, homogeneous models addressed in this pa- 
per. We recall the classical framework and the corresponding loop quantization, including 
a quantization of inverse triad components and general aspects of the Hamiltonian con- 
straint equation. Here we also discuss the two criteria of admissibility of the continuum 
approximation and local stability. 

In section ITTT) we present a quantization of the Hamiltonian constraint. Since the Hamil- 
tonian constraint now also has a potential term depending on the spin connection, a quan- 
tization of the spin connection is required. This is done using the quantization of inverse 
triad components. We derive the (partial) difference equation satisfied by the physical states 
and show that in the continuum approximation, to the leading order, the approximating dif- 
ferential equation is precisely the Wheeler-DeWitt equation with a specific factor ordering 
dictated by the underlying loop quantization. We also show that the requirement of local 
stability is satisfied. We demonstrate that the singularity avoidance mechanism found in 
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the case of Bianchi I |6( continues to hold for these more general models. We conclude that 
the proposed quantization of the Hamiltonian constraint is satisfactory with respect to all 
our requirements. 

In section I1V1 we focus attention on the modification of the potential implied by the 
quantization of the spin connection. The non-trivial and non-perturbative behavior of the 
inverse triad components may be expected to lead to substantial modifications of the po- 
tential at small triad components. In view of the conjectured central role of the Bianchi 
IX potential, we present a brief description of this potential and its implications for the 
modified approach to classical singularities. 



II. DIAGONAL BIANCHI CLASS A MODELS 



In the following, we restrict attention to the so called diagonal, Bianchi class A models. 
Among the class of homogeneous models i.e. models whose symmetry group acts transitively 
on the spatial manifold, is the sub-class of the Bianchi models for which the symmetry 
group contains a subgroup with simply transitive action on the spatial manifold. The simply 
transitive subgroups are classified in terms of three integers, n 1 , parameterizing the structure 
constants as 

Only class A models for which Cjj = admit a canonical formulation. 

In the connection formulation, Bianchi class A models are described in terms of the 
invariant connections and densitized triads given by, 

with constant and canonically conjugate 0}, p{, where u 1 are left-invariant (with respect 
to the subgroup acting simply transitively) 1-forms on the spatial manifold £ and Xj are 
corresponding dual density-weighted vector fields. The 1-forms satisfy the Maurer-Cartan 
equations 

cL/ = -|CVj j Au K . 

In the metric formulation, diagonal homogeneous models are those for which the spatial 
metric can be written in a diagonal form such that only three degrees of freedom remain. In 
the connection formulation the situation is analogous. To begin with one has nine degrees 
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of freedom, 0} but with a three parameter gauge freedom of SU(2) rotations of the index 
i. A further freedom, depending upon the parameters n 1 , would have been available if the 
left-invariant 1-forms were not thought of as a background structure in parameterizing the 
invariant connections. To get the same three degrees of freedom as in the diagonal metric 
formulation, one then restricts the 0} to a "diagonal form", <p) := C(i)A) and correspondingly 
p\ := p T A.\ where the SO(3)-matrix A includes gauge degrees of freedom (see p for details). 
Such a restriction specifies the diagonal models in the connection formulation. This is not 
a symmetry reduction in the same sense as restriction to isotropic models is. Nevertheless, 
these restricted models can be analyzed by following procedures similar to those in the 
context of symmetry reductions. As in reductions to isotropic models, this leads to significant 
simplifications in such diagonalized models because the gauge parameters in A and (almost; 
see below) gauge invariant parameters, cj, are neatly separated 



A. Classical Framework 



The basic variables for diagonal Bianchi Class A models are specified via 

A\ = ctuAjwJ (1) 

in terms of the 'gauge invariant' coefficients cj with the pure gauge degrees of freedom 
contained in the S'0(3)-matrix A. The components cj are not completely gauge invariant 
but subject to residual discrete gauge transformations which change the sign of two of the 
three components simultaneously. A diagonal densitized triad has the form 

Et=p^hiX a j (2) 

where Xj are left-invariant densitized vector fields dual to uj 1 . Being an SO(3)-matrix, A 
satisfies 

A/A/ = 5j , e ijk XjA j jA k K = e IJK 

The triad components p 1 are subject to the same residual gauge transformations as the 
connection components cj to which they are conjugate with basic Poisson bracket 

{ C/ ,/} = 7 «# (3) 

where 7 is the Barbero-Immirzi parameter and k = 8nG the gravitational constant. We 
will use the value 7 = lc ^9 ~ 0.13 fixed by the black hole entropy calculations 
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A diagonal co-triad has the form e\ = a^AjU^ with 

P 1 = \a 2 a 3 \ sgn(ai) , p 2 = \aia 3 \ sgn(a 2 ) , p 3 = \a x a 2 \ sgn(a 3 ) . (4) 

Note that the components p 1 as well as aj can take negative values, but only the overall sign 
sgn(p 1 p 2 p 3 ) (i.e., the orientation) and the absolute values |j/ 1 are gauge invariant. With this 
form of the co-triad, we obtain, in fact, a diagonal metric 

ds 2 = etfju'w' = Y,cli(lu 1 ) 2 . 

i 

The extrinsic curvature also has diagonal form with components Kj = which appear 
in the relation cj = Tj—jKj defining the Ashtekar connection components cj in terms of the 
spin connection components Tj and the extrinsic curvature. The form of the spin connection 
introduces a dependence of the framework on the particular model via the structure constants 



r, = l(*V + ^n*--^-n') (5) 
z \a K a j aja K J 

(K J J K \ 

—n J H — —n K — iy nI J ^ or (-^j J> an even permutation of (1, 2, 3). 

Note that in contrast to the full theory, the spin connection is a covariant object within a 
homogeneous model since coordinate transformations have to respect the symmetry. Con- 
sequently, if non-zero, it cannot be made small by choosing appropriate local coordinates. 
The only model which has identically vanishing spin connection is the Bianchi I model with 
n 1 = 0. Otherwise, the spin connection is non-zero and can even depend on the triad (it is 
a constant r = | in the closed isotropic model, a submodel of Bianchi IX, with n 1 = 1 and 
ai = a 2 = a 3 ). 

The structure constants also appear explicitly in the Hamiltonian constraint 

h = 2k" 1 { [( C2 r 3 + c 3 r 2 - r 2 r 3 )(i + r 2 ) - n l Cl - 7 - 2 c 2 c 3 ] ai 

+ [(cir 3 + c 3 T 1 - rir 3 )(l + 7 -2 ) - n 2 c 2 - 7~ 2 cic 3 ] a 2 
+ [( Cl r 2 + c 2 T 1 -T 1 T 2 ){1 + 7 " 2 ) - n 3 c 3 - j~ 2 Cl c 2 ] a 3 } (6) 
= 2k- 1 [(r 2 r 3 - n 1 T 1 )a 1 + (T 1 T 3 - n 2 T 2 )a 2 + (r\r 2 - n 3 T 3 )a 3 

— \(aid 2 a 3 + a 2 did 3 + a 3 aia 2 )] . (7) 
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In order to derive the classical field equations (for which we consider only positive ai and 
p 1 ) it is advantageous to transform to new canonical variables 

ir i := 2Kjp^ = —aja^ a\a 2 a 3 = —(log a/)' and (8) 
q 1 := |logp 7 such that {q 1 , ttj} = nSj (9) 

where the prime denotes a derivative with respect to a new time coordinate r related to t by 
dt = aia 2 a 3 dr (corresponding to the lapse function N = aia 2 a 3 ). With these new variables 
we have {717, p J } = —2np^5j and 

kNH = Ka ia2 a 3 H = 2 ^(rjj - n 3 r 3 ) + P y (r\r 3 - n 2 v 2 ) + P y (r 2 r 3 - n^x) 

-|(7Ti7r 2 + 7Ti7r 3 + 7r 2 7r 3 )] (10) 

(" 1 ) 2 (? ! ) 2 + (" 2 ) 2 (5 ! ) 2 + (" 3 > 2 ($ ! ) : 

-nV(p 3 ) 2 - nVy) 2 - n 2 nV) 2 - K^tt, + + 7r 2 7r 3 ) (11) 



1 U 1\2 4 1 / 2\2 4 1 / 3\2 41 1222 1322 2322 

= ± [(n ) a\ + [n ) a 2 + (n J a 3 J — n n a\a 2 — n n a\a 3 — n n a\a 3 

-|(7Tl7r 2 + 7Ti7T 3 + 7T 2 7T 3 ) . (12) 

Now one can easily derive the equations of motion, 



which yield 



(loga,)" = 


— 71 j ~ 




aia 2 a 3 H} 


i(logai)" = 


/ 2 2 

(n a 2 


- n 3 a 3 




i(loga 2 ) ff = 


(n x a 2 


- n 3 a 3 


f ~ (n 2 T4 


| (log a s ) w = 


(n x a 2 


2 2 

— n a 2 





(13) 



or, using (g 1 )" = -|(tt 2 + tt 3 ) 

(g 7 )" = -4a 2 /rV . 

For the Bianchi I model we have n 1 = so that log aj = a/(r — tqj). This implies the 
Kasner behavior ai oc t a/ where t = e T and the constraint requires = a\a 2 + a\a 3 + a 2 a 3 = 
i((«i + a 2 + a 3 ) 2 — a 2 — a\ — a 2 ). The coefficients aj can be rescaled by choosing a different 
i(r) which can be fixed by requiring the conventional parameterization a\ + a 2 + a 3 = 1 = 
af+al+a^. As usual, these equations can be solved only if one coefficient, say cti, is negative 



and the other two are positive. Correspondingly, one direction, the first, contracts whereas 
the other two expand toward larger time. When we approach the classical singularity at 
t = 0, space shrinks only in two directions while the third one expands unboundedly. The 
total volume however continues to approach zero according to 0102^3 oc t ai+a2+a3 = t. In 
it has been shown that the Kasner singularity disappears when the model is quantized 
along the lines of loop quantum cosmology. For this result it was important that quantum 
geometry is based on densitized triad variables rather than co-triad variables. In those 
variables, all p 1 cx t l ~ ai decrease to zero when one approaches the classical singularity. By 
contrast, for the minisuperspace described in terms of the scale factors or co-triad variables 
(aj), the singularity is reached when at least one of them goes to infinity. We paraphrase 
this by saying that the classical singularity lies in the interior of the triad minisuperspace 
(containing triads of both orientations and also degenerate ones). 

Bianchi models other than Bianchi I have non- vanishing structure constants and thus the 
evolution of the three triad components can be described as motion in a non-trivial potential 
given by 

w(p\p 2 , P 3 ) = 2 {py(r\r 2 - n 3 r 3 ) + P y(r 1 r 3 - n 2 v 2 ) + P V(r 2 r 3 - n 1 v 1 )} 

= (P l ??l + (P 2 ) 2 r1 + (P 3 m - (p 1 ^ + p 2 T 2 + p 3 T 3 f (14) 

which has infinite walls at small p 1 owing to the divergence of the spin connection components 
(see Fig. |2J). The evolution can then be described approximately as a succession of Kasner 
epochs with intermediate reflections at the potential Q]. For the Bianchi IX model the 
reflections never stop and the classical evolution is believed to be chaotic jl^j ]. 



B. Loop Quantization 

Diagonal homogeneous loop quantum cosmology 6] is first formulated in the connection 
representation where an orthonormal basis is given by the p 7 -eigenstates 

\mi,m2,m 3 ) :— \mi) <g) |m 2 ) ® |m 3 ) (15) 

with 

exp(fimc) 

(c\m) = f . (16) 
V2 sm(ic) 
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(The full Hilbert space of loop quantum cosmology is non-separable, but for our purposes it 
is sufficient to use the separable subspace spanned by the states (JTSJ; see for a detailed 
discussion.) The eigenvalues of the triad operators can be read off from 

p I \m 1 , m 2 , m 3 ) = ^l^m^mi, m 2 , m 3 ) . (17) 

Using the basic operators p 1 one can define the volume operator V = \/\p 1 p 2 p s \ which will 
be used later. Its eigenstates are also |mi,m 2 ,m 3 ) with eigenvalues 

V(m 1 ,m 2 ,m 3 ) = (|t4)^ V\ m i m 2 m 3\ ■ (18) 

A kinematical state \s) is described in the triad representation by coefficients s mi , m2)m3 
defined via, 

\s) = ^2 s m 1 ,m 2 ,m 3 \m 1 ,m 2 ,m 3 ) (19) 

mi,m2,m3 

For a state to be gauge invariant under the residual gauge transformations, the coefficients 
s mi ,m 2 ,m 3 have to satisfy 

^mi,m2,m3 mi,— 7712,1113 ^mi,—m,2,— 1713 m\,m%,— mz • (^^) 

These states are left invariant by the gauge invariant triad operators \p I \ and the orientation 
operator sgn^p 2 ^ 3 ). In calculations it is often easier to work with non-gauge invariant 
states in intermediate steps and project to gauge invariant ones in the end. 

Together with the basic derivative operators p 1 we need multiplication operators which 
usually arise from (point) holonomies hj = exp(c(/)A}Tj) = cos(|cj) + 2A}rj sin(|c/) with 
action 

cos(|ci)|mi, m 2 ,m 3 ) = \{\m x + 1, m 2 , m 3 ) + |mi - 1, ra 2 , m 3 )) (21) 
sin(|ci)|mi, m 2 ,m 3 ) = -\i(\mx + 1, m 2 , m 3 ) - |mi - 1, m 2 , m 3 )) (22) 

and correspondingly for c 2 and c 3 . 

Inverse triad operators 

From the basic operators we can build more complicated ones. We will later need a 
quantization of the spin connection which is a composite operator containing several triad 
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operators. In particular, it also contains inverse powers of triad components which classically 
diverge at the singularity. Since the triad operators have a discrete spectrum containing 
zero, they do not have an inverse. However, general methods of quantum geometry and loop 
quantum cosmology Q] imply that there exist well-defined operators quantizing inverse 
triad components. To obtain these operators one makes use of a classical reformulation, e.g. 

IpT 1 = *T 2 *- a ({ci, V^}) 2 = 67~ 2 ^ 2 (tr A\tMK\ V¥\}) 2 
which can then be quantized to 

j^p = _36 7 - Vb'(j + l)(2j + l)]- 2 (tr, A\tMK\ V¥\]) 2 • (23) 

Here we have indicated that there are quantization ambiguities [3] when one quantizes 
composite operators. The one, most relevant for this paper, is indicated by the subscript 
j of the trace and corresponds to the choice of representation while writing holonomies as 
multiplicative operators. 
This operator acts as 

|p / |~ 1 |mi,m 2 ,m 3 ) := (^£l)~ 1 J\f~ 2 f j (mi)\mi,m 2 ,m 3 ) where, (24) 



(m/) := < k\/\mi + 2k\ > and 
[k=-j ) 

_ j(j + l)(2j + l) 
3 - 3 

The discrete values fj{m) decrease toward lower values for m < 2j [19]. Thus, one can see 
that the classical divergence of the inverse of at vanishing p 1 is explicitly absent in the 
quantized operator. This will be seen to have further consequences for the approach to the 
classical singularity. Furthermore, the state |m 1; m 2 , m 3 ) with rrij = (on which the classical 
inverse triad would diverge) is actually annihilated by due to fj(0) = 0. This allows 

us to define the inverse triad operator (not just its absolute value) by {p 1 )] 1 '■= sgn^)^ 7 ]" 1 
without ambiguity in defining the sign of zero. These inverse triad operators will be used 
later to find well-defined operators quantizing the spin connection. Their eigenvalues are 

(p I )j 1 \mi,m 2 ,m 3 ) = (\^^y 1 U^ 2 sgn{m I )f j (m I )\m 1 ,m 2 ,m z ). (25) 
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2. The Hamiltonian constraint 



Since the classical Hamiltonian constraint contains the structure constants, it must have 
a mode, dependent q uantization. The procedure in the full theory Q motives to ,„ode, 
the holonomy of a closed loop by using the combination hfhjhj 1 ^ 1 in order to quantize 
connection components in the constraint ^J. This works in fact for the Bianchi I model, but 
it must be modified for models with non-vanishing intrinsic curvature. In these models a 
curve which is formed from four pieces of integral curves of vector fields Xj generating the 
symmetry action, which wou.d have holouomy fc,M7 V, do- not dose to form a loop. 
One needs a correction which in jj] has been done by adding a holonomy for the remaining 
curve multiplicatively. However this does not result in an admissible evolution for the closed 
isotropic model derived from Bianchi IX (an admissible quantization for this model has 
recently been presented in [20]). We thus have to find a quantization which satisfies certain 
criteria for a well-defined evolution. 

To appreciate the criteria we recall that the general form of the Hamiltonian constraint 
yields a difference equation for the wave function in the triad representation because of the 
action of holonomy operators Q]. This fact is directly related to the discreteness of quantum 
geometry. However in the classical regime of large volume and small extrinsic curvature, the 
geometry is the familiar continuum (metric) geometry. Thus the difference equation must 
admit solutions which can be very well approximated by solutions of differential equations, 
namely the Wheeler-DeWitt equation which comes from the (Schrodinger) quantization of 
continuum geometry. Such solutions are slowly varying (in the classical regime) and can be 
interpolated by slowly varying solutions of the Wheeler-DeWitt equation and are termed 
as pre- classical. Since quantum geometry depends on the Barbero-Immirzi parameter while 
the continuum geometry does not, absence of 7 dependence in arriving at the differential 
approximation provides non-trivial constraints on the coefficients of the difference equation. 
Admissibility of a continuum approximation is our first criterion. 

However, admissibility of a continuum approximation alone is not enough. Typically the 
order of the difference equation is larger than that of the Wheeler-DeWitt equation. Conse- 
quently there are solutions of the difference equation which are not interpolable by solutions 
of differential equation. Therefore, if in some local region in the classical regime we have a 
"small" admixture of these extra solutions (small so that the continuum approximation is 
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valid), there is no guarantee that contribution of the extra solutions remains small under 
subsequent evolution, i.e. the continuum approximation may get invalidated. That this does 
not happen is our second criterion. It is thus necessary for the stability of the continuum 
approximation. 



A detailed formulation of these criteria is given in [11 1. As a first step, one divides the 
classical regime (large eigenvalues of the triad operators) into smaller cells (local regions) 
which are still large compared to the order of the difference equation. In each of these 
cells, the coefficients of the difference equations are approximately constant. This leads 
to approximating the difference equation by an equation with constant coefficients. In the 
case of ordinary difference equations that result in the isotropic models, the solutions are 
controlled by the roots of a polynomial. The admissibility of continuum approximation puts 
a set of conditions on these constant coefficients which translate into conditions on the roots. 
For the isotropic case in particular, the criterion requires existence of a root with value 1 
and multiplicity 2. The second criterion of local stability further requires that all roots must 
have absolute value less than or equal to 1 in order to prevent exponential growth of the 
non-pre-classical solutions. A further complex conjugation property of the coefficients of the 
difference equations then implies that all roots in fact must have absolute value 1. 

The quantization of the Hamiltonian constraint presented in the next section will be 
shown to satisfy both these criteria. 



III. HAMILTONIAN CONSTRAINTS FOR BIANCHI CLASS A MODELS 

As noted in the introduction, models with non-vanishing intrinsic curvature must have 
non-zero spin connection which, thanks to homogeneity, cannot be made to vanish. There- 
fore, the Ashtekar connection, a linear combination of the spin connection and the extrinsic 
curvature, is not necessarily small even in a semiclassical regime corresponding to small 
extrinsic curvature. While the Hamiltonian constraint in quantum geometry is formulated 
naturally in terms of the connection, for a semiclassical approximation it is desirable to per- 
form an expansion in the small extrinsic curvature. This can only lead to the correct result 
if the spin connection also appears explicitly in the constraint. We will first propose the 
quantized constraint operator, show how it reproduces the classical expression and heuristi- 
cally argue how it can admit a continuum approximation with local stability. The properties 
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of the proposed operator will be elaborated in detail in the following subsections. 
The proposed Hamiltonian constraint operator is: 

H = Ai{ 1 llK)- l Y,e IJK ^[ {T%{AJ)h J {AJ){h I {AJ)Y\hj{AJ))- 1 

-2(f jfj - n L T L )AjA,j} [h K (A) [h K \A) } V 



UK 



where 



h,(A,T) 



1 



p J n K {p K )- 1 + p K n J {p J )- 1 - p J p K n I {{p 1 )- 1 ) 2 



(26) 

(27) 
(28) 



and the {p 1 )' 1 are defined in equation (|25|). Since the inverse triad operators commute with 
the triad operators, there are no ordering ambiguities in the spin connection ([28)1 . Note that 
the holonomy operators h K (A) = exp(c(x)A^Tj) appearing in the second set of braces do 
not need to contain the spin connection since it commutes with V and would drop out. 

In (J26j) the quantized spin connection appears explicitly both in the diagonal 'potential' 
term and in holonomies hi(A^T). For vanishing Ti the constraint operator reproduces the 
one of the Bianchi-I model |4j, |6J. It is easy to see that the term (TjTj — n L T l)AiAj 
comes from the curvature 2-form of the spin connection V for a Bianchi model. In fact, an 
expression via the curvature components would be more general since it could also be used 
for Kantowski-Sachs models. 

The holonomies hj(A,T) are products of holonomies hj(A) and (/i/(r)) -1 . Since we 
have to choose a quantization of the spin connection components T/, which contain inverse 
triads, and also an ordering of the non-commuting A- and T-holonomies, there are many 
more quantization ambiguities compared to the Bianchi-I model. However, qualitative effects 
can often be seen to be independent of such ambiguities. 

To see that the proposed Hamiltonian operator goes over to the classical expression in the 
classical limit, we can ignore the fact that connection components and the spin connection 
do not commute since this will only lead to effects of the order h which disappear in the 
classical limit. The holonomy product hj(A)hj(T)~ 1 can then be written formally as 

h^A^iY)- 1 



which for small Kj 



cos (i( C/ - Ti)) + 2AjTj sin (|( C/ - T 7 )) 
-7 _1 (c/ - Ti) yields 



hi{A)hi{T)- 



2 1 A\r i K I 
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On states ()15|) with large labels mj, the commutator containing the volume operator turns 
into 

h K [h~K,V} \i^l\k K a K 

so that the classical constraint in the form (J7J) is reproduced. 

With regard to the continuum approximation and the local stability conditions, we have 
to take into account the fact that the connection does not commute with the spin connection 
because even small terms can have significant effects should they lead to instabilities. The 
products of holonomies then do not depend on the differences cj — Tj only, as could be 
expected if they were to commute. While ignoring the non-commutativity we have seen 
that expressions like sin (l(c 7 — T/)) appear which, however, are not well-defined in this 
form since there do not exist Cf-operators in quantum geometry (only holonomies of the 
Ashtekar connection are well-defined operators). Motivated by the form of the holonomies 
in (jZ3), 

h^AJ) = (cos(|c / ) + 2A / sin(|c J )) (cos(±f 7 ) - 2A /S in(§f 7 )) 

= \ (e^ Cl e-^ fl + e-^'e^ 1 - 2zA/(e^ c ' e~^ Tl - e'^e^ 1 )) , 

we define sin and cos operators by 

hj(A, T) =: cos (|( C/ - T 7 )) + 2A /S m (|( Cj - T/)) 

with action 

sin (|(ci - Tx)) \m 1 ,m 2 ,m 3 ) := —\i (e^ lCl e~^ Fl - e~^ lcl e^ Fl j \m 1 , m 2 , m 3 ) (29) 

= -^(e- 1 2 iriim < ) \m l + l,m 2 ,m 3 ) 

_ e |ir 1 (m / )| mi _ 1?m25m3 ^ 

and analogously for cos-operators. This means that, by definition, the spin connection 
in trigonometric expressions of cj — Tj, always acts first. In the final expression we used 
eigenvalues Ti(mi) of the spin connection operators, 

Tj\m 1 ,m 2 ,m 3 ) := T J {m I )\m 1 ,m 2 ,m 3 ) . 

Operators of the form ()29|) however, do not produce the basic difference operator in 
the triad representation since they contain the spin connection eigenvalues as coefficients. 
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Nevertheless, for large mi we have in the triad representation for a state \s) 
(|( Cl - rO) s) = \i (e^ r ^ +1 ' m ^s mi+1) 

V / mi ,m2,t(i3 V 



where 



i.e., 



cxp 



,"12,113 

mi,m2,m3 

_p-|«ri(mi-l,m2,m 3 ) 

c s mi-l,m2,m3 

1- -iimir^m/)/- -S 1 "l 



§:=exp{ip I f I /-f£l)s, (30) 

iz^m/r/(mj)^ s mijm2im3 . (31) 

For large mi, therefore, the operators ()29j) do act as basic difference operators on the wave 
function s up to phase factors, and analogously operators cos (|(cj — T/)) act as basic mean 
operators. 

Using the wave function s instead of s, the constraint operator (pjj) contains, besides 
diagonal (potential) terms, a sum of products of a basic difference operator and a basic mean 
operator for two of the three discrete parameters mj. The resulting difference equation for 
s mi ,m 2 ,ra z will thus be of the same form as the difference equation of the Bianchi-I model. 
Consequently, the quantization (j2B| can be expected to yield a locally stable evolution 
equation as in the case of Bianchi-I. We will discuss this in detail in the following subsections. 



A. Discrete evolution equation 



We now return to the constraint operator (|2T)j) and obtain the difference equation in 
detail. From the definitions it follows that, 



hi{A,T) 



% 1 T-l 



g2 1 c~ 2 "h 6 2 6 2 6 2 6 2 — 6 2 62 

+ 2Ai 



Ci - 2iAiSi 



hi(A,r 



-1 6 2 6 2 J -I- £ 2 1 g 2 1 



+ 2Aj 



2/ 



e2 j e 2 J — e 2 6 2 



2/ 



(7] - 2iA 7 5j 



(32) 



(33) 



V" — cos(|cx)l / cos(|c^) — sin(|cx)V r sin(icK) 
-2A# |sin(|c^)\/ cos(|ck) - cos(^c K )V sin(|cx)| (34) 
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Note that the connection components ci do not have a 'hat' on them as these do not 
exist as operators. Their exponentials however are well defined, multiplicative operators. 
For notational convenience these 'hats' have been suppressed. The 'sin', S, has been defined 
without the factor of i also for future notational convenience. 

Then the trace in (|26|) gives (suppressing the hats on the operators): 

27~ 2 {CiCjC\c\ + SjCjSjcl + CjSjCjsl - SjSjSjsl] 

x jy — cos(^ck)V cos(^ck) — sin(~c#)Fsin(|cRr) j 

-2 7 ~ 2 e IJK [SjSjClcl - dSjS\&j + SjCjClS^j + ddS\s\} 
x ^sin(±c K )V cos(±c K ) - cos(~c K )V sin(~cx) | 

-gijk {T/rj - n K T K ) {sm(\c K )V cos{\c K ) - cos{\c K )V sin(\c K )} (35) 
Doing the sum over UK finally yields, 

( ^H tt j c n c^r 1 ^ _i_ Q c ct/^t j_ c /^t ct c c ct ct 
— — — \ n — < L/1O2O1O2 t OiC-2>5i(_/2 1 O1O2O1O2 — O1O2O1O2 
8t J I 

>^ /~i q (~< ct/^t r< c rA ct 1 q c ct ct l 

x jV" - cos(fc 3 )Vcos(|c 3 ) - sin(fc 3 )V sin(|c 3 )J 
I c o rA n q ct/^t \ q n cA eft a_ n n ct ct 

— S Oi/J2>^x^ / 2 — , - y l | J2 | Jl t - / 2 | Jl t - / 2 , - y l >- , 2 ' ( - / l ( - / 2*-'l ^2 

_i_ c c rA rA fi a q\ rA \ q n rA oA \ H r ct ct I 

"T>-'2>-'l < - / 2^1 — 0201020-^ 020iLy 2 >-) 1 -f- C^Oi^Oj f 

x |sin(|c3)V"cos(ic 3 ) - cos(§c 3 )V sin(|c 3 )| 

-7 2 {rir 3 - n 3 r 3 } jsin(fc 3 )V cos(|c 3 ) - cos(~c 3 )Vsin(|c 3 )j 

+ cyclic (36) 

The terms containing V and f / are both diagonal when acting on the basis vectors 
\mi,m 2 ,m 3 ). They give: 

T I \m 1 ,m 2 ,m 3 ) = r I (m 1 ,m 2 ,m 3 )\mi,m 2 ,m 3 ) where, (37) 
T I (m 1 ,m 2 , m 3 ) := ^A/"" 2 (mjn K fj(m K ) + m K n J fj(mj) - mjm K n l ' Nf A ! /|(m 7 )) 

using the notation of (J21J) . 

Next, 

(v - cos(|c 3 )y cos(±c 3 ) - sm(§c 3 )Vsin(§c 3 )) \m 1 , m 2 , m 3 ) = 

V(mi, m 2 , m 3 ) - ~V(mi, m 2 , m 3 + 1) - ~V(mi,m 2 , m 3 - 1)^ \mi,m 2 , m 3 ) (38) 
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(sin(±c 3 )Vcos(±c 3 ) -cos(§c 3 )Vsin(ic 3 )) \m 1 ,m 2 ,m 3 ) = 

- (V(m l ,m 2 , m 3 + 1) - V(m 1 , m 2 , m 3 - 1)) |mi, m 2 , m 3 ) (39) 
The action of the Cj, 5/, C 1 ), Sj on the basis states can be expressed as: 

Ci|mi,m 2 ,m 3 ) = ^ e^ 61 ^™ 1 '™ 2 '™ 3 Vi + d, ™2, ™3> (40) 
2 ei =±i 

5117711, 7712, 7713) = \ ^ e x e -*«i r i("»i."*.™3 >|mi + 6i, m 2 , m 3 ) (41) 

2 £1 =±1 

Cjlmi,^, ™3> = \ e-^ Vl{mi+ ^ m ^ m3) \m 1 + e' l ,m 2 ,m 3 ) (42) 

e[=±l 

S\\ mi ,m 2 ,m 3 ) = ~ e 'i e-y^ {m ^ m ^\ mi + e'^m^m^ (43) 

6' 1= ±1 

In each group of four terms, the phase factor will be the same since each consists of two 
C, S type of terms and two of C\ types of terms. We will just have different factors of 
the 61,61- Thus, for example, 

{CiC 2 CM + S&Slcl + C^Clsl - S 1 S 2 S\S^ |mi, m 2 , m 3 ) = 

\ \ r g-f {(ei+e' 1 )ri(mi+e' 1 ,m.2+e2> ,n 3) + £2^(7711 +e' 1 ,m2+e 2 !" J 3) + e 2 r 2( m i: m 2+e2> m 3)} 

16 ^ 

ei,e2,e'i> £ 2 

x (1 - eie'i - e 2 e' 2 - eie / 1 e 2 e / 2 ) |toi + ei + ei, m 2 + e 2 + e' 2 , m 3 ) , (44) 
\^C 2 C\C\c\ + iS' 2 CijS' 2 Cj + C 2 jS'iC2jS'j — iS' 2 jS'ijS'2jS'J ^ |mi,m 2 ,m 3 ) = 

J_ \ r g-f {(e2+e2) r 2( m l+ ,E 'i. m 2+<:2' ,n 3) + eilMmi+e'^rr^+e'^m^) + e' 1 r l( ra l+ e 'i) m 2,™3)} 

16 ^ 

ei,e2,e'i> e 2 

x (1 - eie'i - e 2 e' 2 - eie / 1 e 2 e / 2 ) |mi + ei + e'i, m 2 + e 2 + e' 2 , m 3 ) , (45) 
\S X S 2 C\C\ - C X S 2 S\C\ + S x Cid[S\ + CiC 2 SM} |m 1; m 2 , m 3 ) = 

J_ \ r e -3{(ei+e'i) r l( m l+ e 'l> m 2+«2>"i3) + ^F 2 {m\ +e' 1 ,m2+e' 2 ,m 3 ) + e 2 r 2 (mi,m 2 +e 2 ,m 3 )} 

16 ^ 

ei,e 2 ,ei,e 2 

x (e x e 2 + e / 1 e 2 - eie' 2 + e[e 2 ) \mi + ei + e\, m 2 + e 2 + e' 2 , m 3 ) (46) 

{StStClCl - C^Slcl + S^Clsl + CaCiS'JS'l} |m 1? m 2 , m 3 ) = 

18 



and 



_ \ r „-f{(<:2+e2) r 2(mi+e' 1 ,m2+e2: m 3) + eil 1 ! (7711+6^,7712 +4 > m 3) + ejri(mi+ej,m2,m3)} 
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ei,£2,e' 1 ,e 2 

x (eie 2 - + eie' 2 + eie 2 ) |mi + ei + ei, m 2 + e 2 + e' 2 , m 3 ) . (47) 
It is convenient to define different basis vectors at this stage, namely, 

|mi,m 2 ,m 3 ) := e^' miri{mi ' m ^ m3) \m 1 '^^m 3 ) (48) 

corresponding to the transformation (|30p. This introduces additional phases which will turn 
out to be useful in our discussion of the continuum approximation. The transformation to 
the triad representation via, 



\s) := 22 S rn 1 ,m 2 ,m 3 \m 1 ,m 2 ,m 3 ), (49) 
777,1,7712,7713 

leads to an equation for the s miirn2jm3 : 

61,62,4,62 

+ ^2 7~ 2 ^23(^i,^2,m 3 ;e2,e3,e 2 ,e 3 )s miim2 _ £2 _^ im 3_ e 3_^ 

62,63,e 2 ' e 3 

+ ^2 7~ 2 ^3l(^l,m2,m3;e3,ei,e 3) e/ l)5m 1 - ei -6' 1 ,m2,m3-63-6! i 
63,61,63,61 

+ B(m 1 ,m 2 ,m 3 )s miim2jm3 , (50) 
where the coefficients are rather complicated expressions given by: 

vf 2 (m 1 ,m 2 , m 3 ; ei, e 2 , e[, e' 2 ) := V{m x -e x - e[, m 2 - e 2 - e' 2 , m 3 + 1) 

± V(mi - ex- ei, ma - e 2 - e' 2 , m 3 - 1) , (51) 

«i 2 (mi, m 2 , m 3 ; e x , e 2 , ei, e' 2 ) := V{m x -e x - ei, m 2 - e 2 - e' 2 , m 3 ) 

- ^vf 2 (mx,m 2 , m 3 ; ei, e 2 , ei, e 2 ) , (52) 



-'12,1 



(mi,m 2 ,m 3 ;ei,e 2 ,ei,e 2 ) := - ^ m I r I (mx,m 2 , m 3 ) 

+ (ei + ei)ri(mi - ei,m 2 - e 2 ,m 3 ) 
+ (e 2 + e' 2 )r 2 (mi - ex, m 2 - e 2 , m 3 ) 
+ ei { Tx{mx - ei,m 2 - e 2 - e' 2 ,m 3 ) 

- T^rni - ex,m 2 - e 2 ,m 3 ) } (53) 
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and 



(j)i 2 , 2 (m 1 ,m 2 ,m 3 ;e 1 ,e 2 ,e' 1 ,e 2 ) := - ^ m I T I {m u m 2 , m 3 ) 

i 

+ (ei + ei)ri(mi - e 1 ,m 2 - e 2 ,m 3 ) 
+ (e 2 + e 2 )r 2 (mi -e 1 ,m 2 - e 2 , m 3 ) 
+ e' 2 { r 2 (mi - ei - ei, m 2 - e 2 , m 3 ) 

- r 2 (mi - ei,m 2 - e 2 ,m 3 ) } (54) 

With obvious abbreviations, e.g., (m 1 ,m 2 ,m 3 ) = rh, (ei,e 2 ,e f 1 ,e' 2 ) = e 12 and 
A 12 (m 1 ,m 2 ,m 3 ; ei,e 2 ,ei,e 2 ) = A 12 (m; e i2 ) etc. , we have 

Ai 2 (m; e 12 ) := m 12 (m;e 12 ) (1 - - e 2 e' 2 - e^e.e,) ^-^(to) - e -§<^,i(™A 



and 



+ ^r 2 (m; ei 2 ) {e-i^^^) ( ei e 2 + e'^ - e 1 e' 2 + eie 2 ) 

+ e -^».i(«;ei2) ( ei62 _ e ; e2 + eie ' 2 + e ; e ' 2 )| 

£(m) := 8 (e^ m ' r '<™)) 



(55) 



x 



(56) 



v 12 (m; 0) {ri(m)r 2 (m) - n 3 r 3 (m)} 
v 23 (m;6) {T 2 (m)T 3 (m) - n^m)} 
v^fh-jO) |r 3 (m)ri(m) — n 2 r 2 (m)} 

— ?K^f 2 ( ' P %'El m l r rt™)\ C/matter 

The matter Hamiltonian acting on the matter field dependence of the wave function has 
also been incorporated in the coefficient B{m). 

Noting that the coefficients ui 2 ,Ai 2 ,B involve vf 2 ,V and that the volume eigenvalues 
have a separable product structure, we can absorb the volume factors in the wave functions 
as follows. Define, 



m 1 ± y/1 — m~ 



if m 7^ 
if m = 



Then, 



vf 2 (m; e) = V(m 1 — e± — e[, m 2 - e 2 - e' 2 , m 3 )rf ± (m 3 ) 

, _ T , / \ f , <^ + (^3) ] 

u 12 {m;e) = V{m\ — e\ — e 1 , m 2 — e 2 — e 2 , m 3 ) < 1 > 



(57) 



(58) 
(59) 
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Now defining i muma>ma : = V(m)s mumatma and p m ^ r 2 , m3 := H%^ >m3 /V(m) for 
V(mi,m2,ms) 7^ 0, one sees that all volume eigenvalue factors are absorbed in the wave 
function. This preserves the gauge invariance condition on the S but now in addition, due 
to the explicit volume eigenvalues, we also have that t mi>m2jm3 = if any of the m^'s equal 
zero. 

For Bianchi I, the spin connection is zero. The phases 0i2,i, 0i2,2 ; • • • vanish. The coeffi- 
cient B reduces to the matter term only, while A\ 2 reduces to : 

A 12 (rh; e u ) = d~(m 3 ) (e^ + e' x e' 2 ) (60) 

Equation (jSOjl then becomes, 

d (7Ttl) (tmi ,m2+2,m3+2 tm\ ,rri2 — 2,7713 — 2 ^7711,7712+2,7713— 2 ^7711,7712 — 2,7713+2) 
-\- d (7712) (^7711+2,7712,7713+2 "I - ^"m\ — 2,7112,7773 — 2 ^mi- 2,7712,7713+2 ^7771+2,7712,7773 — 2) 
-|- (^3) (^mi+2,m2+2,m3 "I - ^mi- 2,7712— 2,m3 ^7711+2,7712—2,7713 ^7711—2,7712+2,7713) 

,7712,7713 ^7711,7712,7713 

(61) 

which matches with the equation derived in 



B. Continuum approximation 

We will first derive a simplified difference equation which approximates the exact one for 
large m;. Then using an interpolating function, we will show that the approximate difference 
equation gives a differential equation for the interpolating function, which matches with the 
Wheeler-DeWitt equation at the level of leading terms. This will verify that the continuum 
approximation is admissible. 

A/"? 

Consider first mi ^> 1. Then, for k <€. m, fj(m + k) — > - 1 - which is independent of k to 
the leading order. This implies that 

r 1 / K m J . j m K rmjiriKl 

Tj ~ -<n Yn J n =— > . (62) 

2 ^ ttik trij trtj J 

The phases </>i2,i, ^12,2 etc. all become equal to the common phase appearing in the coef- 
ficient B and cancel out. The coefficients Ajj become exactly the same as in Bianchi-I with 
d~(mj) rs mj 1 and the equation becomes, 

ei,e2,£'i,e 2 
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+ 8 m 3 1 (FiTa - n 3 T 3 ) t mi , m2 , m3 + cyclic (63) 

To arrive at a differential equation, we look for a function T(p 1 ,p 2 ,p 3 ) such that 
tm 1 ,m 2 ,m 3 '■= T(p 1 (m\ ) , p 2 (7712) , p 3 (m.3) ) satisfies the approximate difference equation. Here, 
:= T^^m/ and thus £p 7 = is the change induced when the m/ change by 1. In 
particular, this implies: 

i mi+kum2+k2!m3+ka := f(p I (m I + k I )) = f(p I (m I ) + k I Sp I ) 

fyf i <9 2 T 

- r (pV/)) + E + 5 £ 

+ ■ • ■ (64) 

It follows that terms of order (dp 1 ) , (dp 1 ) 1 cancel out leaving a second order differential 
expression to the leading order. We also have the pre-factors •y~ 2 mj 1 which can be eliminated 
using mi = %-7 _1 . The spin connection then becomes just the classical expression without 
any factors of 7^ 2 . All factors of 7 cancel out. Multiplying by ^p 1 p 2 p 3 the final equation 
becomes: 



1 2 d 2 T(p 1 , p 2 , p 3 ) 1 2 3 ~ 12 3 

~l P P dp 1 dp 2 P P ^ ~ n ' ^ ,P ,P 



cyclic 



= _l Kp i^0 r ^(piyy)f( p i iI ?y> j (65) 

We have thus verified that the constraint operator (|26|). admits a continuum approxima- 
tion. 

Note that using an interpolating function, T(p l ,p 2 ,p 3 ), one can always arrive at some 
differential expression. That all dependence on 7 disappears from the differential equation 
is a non-trivial property of the locally approximated difference equation and is essential 
for matching with the Wheeler-DeWitt equation which knows nothing about the Barbero- 
Immirzi parameter. Had the phases not canceled, we would have acquired terms from the 
Taylor expansion of the 17, which could have produced lower order differentials with a 7- 
dependence left over, destroying the hope of matching with the Wheeler-DeWitt equation. 
That the spin connection controls the phases in a cancelable manner is thus also a non- 
trivial property of the constraint operator. Finally, the fact that we get a purely second 
order differential expressions is determined by the structure of the coefficients in the locally 
approximated difference equation jll]. This is directly responsible for getting the Wheeler- 
DeWitt equation in a particular factor ordering appearing in (J55j) . 
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The wave function T has been obtained after some redefinitions. Since V(m) — > ^\p 1 p 2 p 3 \ 
in the pre-classical limit, T is related to S as T ~ ^/\p 1 p 2 p 3 \ S. Thus, the Wheeler-DeWitt 
equation for our original pre-classical wave function S is 

" 4 i 2 9VbWI s{p\p 2 ,p 3 ) 



p 

-p p 



+ P y (r x r 2 - n 3 r 3 ) v / bWl £(p\ /A p 3 ) 



+ cyclic 



4 dp 1 dp 2 

K \ p ± p y 1 1 patter (p 1 , p 2 , p 3 ) ^(p 1 , p 2 , p 3 ) . ( 66 ) 

Note that the ordering of an isotropic sub-model is not determined uniquely by using p 1 = 
p 2 = p 3 =: p because p 1 and d/dp 2 , say, commute in the anisotropic model but not after 
using isotropy. However, there is also a unique ordering for isotropic models J^, I2]] which 
can be derived after doing the reduction in loop quantum cosmology 



C. Local Stability 

Our next requirement is about local stability of the evolution equation around pre-classical 
solutions. This is most directly formulated and illustrated in the context of an ordinary differ- 
ence equation ^j] . The equation we have however is a partial difference equation. Although, 
we do get a partial difference equation with constant coefficients in a local approximation, 
their solutions are not as easily analyzable in terms of roots of some polynomial. Neverthe- 
less, the property of local stability can still be analyzed. 

Since we have three m/s, even if all are taken to be large, it is possible to have a subset 
of them to be much larger compared to the remaining ones. To begin with, we consider a 
neighborhood of a point where all mj are large and of the same order; other regions can be 
explored similarly. For 7 we assume the value derived from black hole entropy calculations 
which is of order 10 _1 . Since all m/s are large, it follows exactly as before, that Tj are given 
by the equation (|62|) and all phases cancel out giving the locally approximated difference 
equation ([63)1 . Furthermore, since we assume the mi to be about the same size, the Tj are 
approximately equal and taking the Bianchi IX model for definiteness, « ~. We obtain 

d (?Tli) {t mi ,7712+2, m.3+2 tmi,m2 — 2,m,3 — 2 ^mi,m2+2,rri3 — 2 ^mi,m2-2,m3+2) 
d (vTl2) (^mi+2,m2,rri3+2 tmi — 2,m2,m3~2 ^mi-2,m2,m3+2 tmi-\-2,m2,m^ — 2) 
-|- d (7713) (^mi+2,m2+2,rri3 ^mi-2,m2-2,nj3 ^mi-2,m2+2,nj3 ^mi+2, rri2 — 2,1x13) 

= -^t\p™^ mo t mmmz + 2 7 2 (GT(m 1 ) + d-(m 2 ) + GT(m 3 ))t mi , m2 , m3 . (67) 
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In this equation, mi, m 2 and m 3 all play the same role, but without loss of generality we 
can select m 3 to be the evolution parameter in terms of which we analyze local stability. To 
express it clearly, let us introduce the notation 

r ^i,i^2 — f 1/1/ n -1-9 

which allows us to write the equation as 

d-(mi)(z°J 3+2 - z°J 3 _ 2 ) - d"(mi)(^+ 2 - z^ 2 2 ) 

i A- ( \t 2.° 2,0 \ i-i \( -2,0 -2,0 \ 

+ d {m 2 )(z m3+2 - z m3 _ 2 ) - d {m 2 )(z m3 ' +2 - z ms '_ 2 ) 

j- 1 \l 2,2 , -2,2 2,-2 — 2,— 2\ 

= -d {m 3 )(z m3 + z mi ' -z m3 -z ma > ) 

- 2 7 2 ( 7 P + d~(mi) + d~(m 2 ) + d'(m 3 ))z° m ° 3 (68) 

where 

p ■= K 4p mattcr 

is a small constant (in the local approximation) characterizing the size of the matter contri- 
bution. 

For m/s of the same order m, the equation can be further simplified as 

/ 0,2 0,2 \ _ / 0,-2 _ 0,-2 \ , / 2,0 2,0 \ _ / -2,0 _ -2,0 \ 

\ Z m 3 +2 Z m 3 -2J \ Z m 3 +2 Z m s -2J ' \ Z m s +2 Z m 3 -2J \ Z m 3 +2 Z m 3 -21 

- -(« + z£* - 4 3 " 2 - z ^~ 2 ) - 2 7 2 (3 + imP)zZ (69) 

We are interested in checking if perturbation of slowly varying (along all three directions) 
solutions, continue to remain as perturbations under the m 3 -evolution. To this end, we 
assume that we have some fixed solution which is slowly varying. Then the first bracket on 
the right hand side is small while the second one is suppressed by 7 2 ~ 1CT 2 (mP is small 
even for large m because the energy density decreases at least as m~ 3 ). 

This is a linear equation among four combinations analogous to the equation of a plane 
in R 4 passing close to the origin. It can be solved approximately by introducing three, fixed 
orthonormal vectors, S a , a = 1,2, 3, each of which is orthogonal to the vector (1, —1,1, — 1) 
in 1R 4 . Explicitly, 

3 

(^3+2-^3-2) « E^( m 3)(£ a ) ,2 
a=l 
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3 

(^3+2-^7-2) ~ ^2 Qa(m 3 )(£ a ) ^ 2 
a=l 
3 

(^3+2-^3-2) ~ ^2Qa(m 3 )(£ a ) 2t0 

a=l 

3 

(^3+2 - ^3-2) « E ^( m 3)(^ a )-2,0 (70) 
a=l 

In the above, corrections of the order z^" 2 — Zm^ 2 as well as of the order 7 2 (3 + 7mP)z^ 
are ignored and Q a {tTfi) are three functions of 777,3. The Q a (^3) can be assumed to be small 
for z to represent a solution close to a pre-classical one. 

Thus we obtain ordinary, non-homogeneous difference equations with constant coefficients 
for the four z^ 2 . All the equations have the same structure, and their solutions consist 
of a combination of the particular solution and linear combinations of solutions of the cor- 
responding homogeneous equations. The difference between a perturbed solution and the 
fixed solution of course satisfies the homogeneous equations. 

The local stability now requires that these differences do not grow exponentially. Equiva- 
lently, the characteristic roots for the homogeneous equations must have absolute value less 
than or equal to one. 

In our case, the characteristic polynomials are the same for all four equations namely, 
z A — l =0, with roots equal to ±1, ±z. These obviously satisfy the condition of local stability. 

Because of local stability, perturbations of slowly varying solutions continue to remain 
close to them under evolution along the m 3 direction. This also ensures self- consistency 
when we consider evolutions along the mi, 777.2 directions in a similar manner by keeping the 
corresponding right hand sides close to zero. 

Had we not neglected the terms on the right hand side of eq. 1)69)1 . we would have got 
an extra, common non- homogeneous term which is 1/2 of the right hand side of eq. ljHH)) . 
Clearly, this does not affect the behavior of the differences. 

To summarize, although we got a partial difference equation, we were able to implement 
the idea behind the local stability criterion by obtaining a system of ordinary difference equa- 
ions. That these equations are non-homogeneous, does not affect the stability properties 

m 
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D. Absence of Singularities 



We have now verified that the proposed quantization (|26|) has a valid continuum approx- 
imation and is locally stable around pre-classical solutions. Thus, it is possible to start 
with initial data at large volume which are close to a continuum solution and to follow 
the evolution to smaller volume, toward the classical singularity. In this regime, quantum 
effects of the discreteness and also ordering issues become more and more important and 
the full equation has to be considered without approximations. Fortunately, it is possible to 
see essential features of the absence of singularities quite generally in a way similar to the 
isotropic case jjj. Since the constraint equation is used as a recurrence relation for the wave 
function, the term with lowest order in the labels mj should have non- vanishing coefficient, 
except when one of the labels vanishes such that the corresponding component of the wave 
function, i.e. its value at singular surfaces in minisuperspace, will drop out (see p for de- 
tails). It can readily be verified that this general behavior is true for the evolution equation 
of an arbitrary Bianchi class A model. The difference equation resulting from (l26l) is more 

n n 

complicated than in the isotropic [SJ or Bianchi I |6] case; but since the exponentiated spin 
connection never vanishes and the additional potential term containing the spin connec- 
tion is diagonal, only the commutator of holonomies with the volume operator determines 
whether or not a lowest order coefficient vanishes. The structure resulting from this term is 
the same as in the Bianchi I case, and thus the evolution is singularity-free in the same way. 

This shows that the wave function can be extended uniquely from the region of min- 
isuperspace with positive orientation to the region with negative orientation and that the 
orientation again provides us with a branch 'before' the classical singularity. 

For Bianchi models other than Bianchi I, however, the classical approach to the singularity 
is more complicated than the simpler Kasner behavior. This is most dramatic in the Bianchi 
IX model where different Kasner epochs succeed each other infinitely many times. It turns 
out, as we will discuss in the next section, that the classical approach to singularities is 
modified at small scales by quantum geometry effects which is expected to lead to a simpler 
picture. 
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E. Quantization Choices 

Even after imposing the condition of local stability the constraint expression is far from 
being unique. There may be other choices for the kinetic term which lead to a locally stable 
behavior, but they would have to be very special and difficult to find. The advantage of the 
quantization presented here is that it is rather simple and provides a general formula for all 
diagonal Bianchi class A models at once without having to find an expression in each case 
separately. 

As for factor ordering ambiguities, the main issue is that of ordering the commutator 
to the right or to the left. This ordering is the same here as in the Bianchi I case and 
also in the full theory. Thus, the special issues one has to deal with in models with non- 
vanishing intrinsic curvature do not change this result. Since the non-singular evolution 
depends significantly on this kind of ordering, it is robust. The main purpose of this paper 
is to prove the existence of one quantization for all models in our class which has the 
correct semiclassical limit, results in a locally stable evolution and is also non-singular in 
the quantum regime. That such a general quantization fulfilling all the different requirements 
exists is a non-trivial fact and provides a further check of the methods employed in the full 
theory. 

The second issue, besides the singularity-free evolution, which is robust under the choices 
we have is the semiclassical behavior and corrections to it. The main effects in this context 
are insensitive to factor ordering and other choices. As we will discuss in the next section, 
there are non-perturbative effects coming from the potential term which lead to drastic 
differences to the classical behavior. Since they come from the potential term and are based 
on general features of quantizations of inverse triad components, they have to be present in 
any quantization. 

IV. MODIFIED BIANCHI IX BEHAVIOR 

The quantization (|26|) of the Hamiltonian constraint yields consistent quantizations of the 
dynamics of Bianchi class A models which in the limit of large triad components reproduce 
the correct classical behavior. At finite values of p 1 , however, there are always correction 
terms to the classical equations of motion (fT3|) which can change the classical evolution. It 



27 



is possible to derive different types of corrections and to include them in effective classical 
equations of motion. There are the following sources of correction terms: First, we can 
keep more terms beyond the leading order in a Taylor expansion of difference operators 
which results in a higher order Wheeler-DeWitt equation. These correction terms would 
correspond to higher curvature terms in an effective action. The second source is the non- 
commutativity of the A-holonomy hj(A) and the T-holonomy hj(T) in (|2*U|) . Since we have, 
for example, {ci, Ti} = / ~fKp 2 p 3 (p l )~ 3 n l which in a quantization would appear as the diagonal 
operator 

zftAi{Sj7} = —hi(A) [/ii(A) _1 ,fi] , 

the coefficients of the Wheeler-DeWitt equation receive corrections of the order o^^/p 1 ). 
Thirdly, we have to use well-defined inverse metric operators when we quantize the spin 
connection components in (|2fi|) which leads to a modified behavior at small scales (as in the 



isotropic case where it implies inflation 22j). 

Since the last effect is non-perturbative while the first two are perturbative (in 7) , it can 
be expected to have the most dramatic consequences. Therefore, we will focus on this type 



22], 



of correction in our discussion of a modified approach to classical singularities. As in 
we exploit the fact that there are quantization ambiguities j3] which can extend the range 
of a modified classical behavior If we choose the ambiguity parameter large enough 

it is reasonable to use effective classical equations of motion to describe the trajectory of 
a quantum wave packet while still experiencing the modified behavior at small scales. In 
particular, the spin connection components will always decrease when we approach the 
classical singularity close enough whereas they always have some diverging components in 
the purely classical scenario. The diverging components lead to infinitely high walls in the 
potential (|14jl on minisuperspace and to an infinite number of reflections in the Bianchi 
IX model leading to chaotic behavior |l5[. With the modified spin connection, the walls 
do not grow arbitrarily high but instead shrink once a small triad component is reached. 
Consequently, one expects that the reflections will cease after a finite number of times and 
the wave packet will simply move through the classical singularity in a final Kasner epoch. 

Here, we will make this explicit by deriving the potential and its changes due to the 
modified spin connection. A detailed analysis of the modified approach to the singularity 



will be presented elsewhere 



24, 



25|. 



The quantized spin connection operators Tj have been defined in equation ()28|) with the 

28 
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FIG. 1: The function F(q) (solid line) compared to its large-g limit q 1 (dashed). 



inverse triad operators denned in equation (|24jl . The functions fj(m) can be approximated 
(with an accuracy increasing with j) in a manner similar to that discussed in [l9(. In terms 



of q :— Yp one obtains 



/i(2j?) 
F(q) 
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fF(q) 



where, 
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~ q^ 1 for g ^> 1 
~ 4q 2 for g < 1 



\q — 1 1 2 



5<7 



(g + l)5 -sgn(g-l)|g-l|5 | (71) 



(see Fig. Q). Limiting behaviors for the eigenvalues of (p 7 ) 7 - 1 are 



(P^/W « (7^7') F(m I /2j)sgn(m I ) 



p I {m J ) 1 



4 (7^pi) 3 p / (m J ) 2 sgn(m 7 -) 



for |p 7 (mj)| 3> 7^|j 
for (^(mj)! <§C 7^pj 



(72) 
(73) 
(74) 



Here j/(mj) = ^•yip'mi are the eigenvalues of p 7 . The eigenvalues of the spin connection 
operators are obtained as, 



Tj (in) « — [ n J rnxFimj sgn(mj) + n K mjF{rriK /^j) sgn(m^) 
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— n I mjmKF 2 (mi /2j)/2j 



(75) 



n 



1 
2 

1 

4? 



j P K (m K ) ^ K p J (mj) 
p J {mj) p K {m K ) 



— n 



jP J (mj)p K (m K ) 
(ptyirm) 



(Jfe » l) (76) 



n J mK^j sgn(mj) + n K mjm z K sgn(mx) — n 



jmjmKm] 



« 2 ( 7 ^rV(™j)p X (™*) [ nVNI+nVWI ] (^7 « l) (77) 

In the last equation we have dropped the sub-leading third term. Evidently, the eigen- 
values of the spin connection operators approach their classical expressions for large values 
of the triad components while for small values the quantum modification shows up. To find 
the potential in this regime, we must use the first line of eqn. ()14|) since the second follows 
by using the classical expression for Tj. For |p 7 | <C •yipj this leads to, 



Wj(p) := 2pY (T^T® - n 3 Ti j) ] 



cyclic 



« [8(7#)-VpV)VVI +n 3 \p 3 \)(n 1 \p 1 \ + n 3 |p 3 |)) 
-4(7^j)~ 3 (pV) 2w3 ( w Vl +n 2 \p 2 \)] + cyclic. 

Dropping the o(j~ 6 ) term compared to the o(j~ 3 ), for |p 7 | <C 7-O', we get 



^3 1 3 



By contrast, for |p 7 | 3> 7^', the potential takes the form, 
Wj(p) « 7 



{pV^j-V } 2 + {pVbV 1 ^ } 2 + {p 3 pV)-V } 2 

- 2 {(p^nV + (p 2 ) 2 nV + (p 3 )Vn 2 } ] 



Specializing to Bianchi-IX this is: 



W 3 {p) 



,_|_ [ (p3) 2 ( | p l| 3 + b 2|3 ) + (p l )2 ( b 2|3 + b 3|3 ) + 

(p 2 ) 2 ( |p 3 | 3 + bT)l 



\P'\ 



p 1 p 2 

p3 



+ 



^ 1 



) 2 -2{(p 1 ) 2 + (p 2 ) 2 + (p 3 ) 2 } 



7 < 1 



> 1 



(78) 



^(p) » -4( 7 ^r 3 [ ™w) 2 ( ipT+ ip 2 i 3 ) + ^W) 2 ( bT+ ip 3 

+ nV(p 2 ) 2 ( |p 3 | 3 +|pT ) ] (79) 



(80) 



(81) 
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FIG. 2: Classical potential as a function of x := 2p x and y := 2p 2 /j£p with 2p 3 /j£p = 10. 




FIG. 3: Modified potential with collapsed walls at small p 1 and p 2 (with j = 3). 

A sample snapshot of the modified potential is shown in Figs. El and 0] together with 
the classical potential (Fig. |2J) for comparison. With the modified potential, the classical 



31 



30 30 



FIG. 4: Modified potential with collapsed walls at small p 1 and p 2 (diagonal canyons) with j = 3. 
Also regions for negative p 1 are shown to which a test particle in the potential can move; with the 
original potential only the lower quadrant is allowed. 

Hamiltonian (fTU|) becomes 

KNHfi^q 1 ) = -^(vr 1 7r 2 + 7r 2 7r3 + 7r 3 7r 1 ) + ^(e V ,e 292 ,e 293 ) (82) 

We have thus derived an effective potential modified by the non-perturbative quantum 
effects. Other correction terms, such as higher order contributions to the kinetic part coming 
from the discreteness, are perturbative and less important as remarked earlier. Note that 
classical regime is identified by rrij = p / (-^ £ ) _1 3> I • The quantum modifications however 
are dominant in the small q = p regimes. Choosing larger values of the ambiguity parameter 
j, thus allows us to access the quantum effects even while staying in the classical regime. We 
can thus study the qualitative effects of quantum modifications within a classical framework 
by using the effective potential instead of using a fully quantum description in terms of 
a wave function on minisuperspace. Effective classical equations of motion, in terms of an 
external time parameter, can be derived from the Hamiltonian (|82j) as we did in Section 111 Al 
The corresponding classical trajectories can be thought of as (approximately) describing the 
motion of a quantum wave packet, giving suggestions regarding the full quantum behavior. 

Alternatively, one can also analyze the modified behavior in terms of an internal time 
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FIG. 5: Classical potential at constant volume V\ = V2 = (57^)2, and V3 = (37^)2 

(from left to right). While moving in the triangular valley the universe is pushed toward the 
classical singularity (cc = y = 0) by reflections on the moving left wall. 




FIG. 6: Modified potential with collapsing walls (for j = 3) at constant volume V\ = (^7^p) 2 , 
V 2 = (57*1)5, and V 3 = (3rfj£)* (from left to right). 

namely the volume. The Hamiltonian (|82|) is in fact hyperbolic. The potential now has a 
non-factorizable time dependence and the analysis of motion is more complicated. Figs. 03 
and El show the classical and the modified potential, respectively, at fixed volume \/\p 1 p 2 p 3 \ 
rather than at fixed p 3 as in Figs. El and E3 One can already expect from the figures that the 
behavior will change dramatically once the universe reaches a small triad component. For a 
detailed analysis see 0, Q • 



V. CONCLUSIONS 



In this paper we have extended the methods of homogeneous loop quantum cosmology 
to models with non-zero spin connection. These models have non-trivial intrinsic curvature 
from the symmetric background which requires a special treatment compared to the Bianchi 
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I model as well as compared to the full theory. Initially, therefore, part of the guidance from 
the full theory is lost and there are more ambiguities when one quantizes the Hamiltonian 
constraint. 

In |h| two conditions which have to be imposed for a reasonable quantization, have 
been introduced and discussed in detail. These stem primarily from the requirement that 
the loop quantization admit physical states which, in a semiclassical regime, have much 
less sensitivity to variations on the Planck scale so that the idealization of the continuum 
geometric formulation is well justified. This has been made more explicit by the formulation 
of the continuum approximation. 

In Section 11111 we introduced a quantization of the Hamiltonian constraint for all diag- 
onalized Bianchi class A models, and showed that the two conditions are satisfied. It was 
important to take the freedom of phases in the triad representation into account which we 
did by transforming to a new wave function s in ([30)1 . The transformation was dictated by 
the requirement of reproducing the Wheeler-DeWitt equation as the leading approximation. 
It is satisfying to note that the same transformation also appeared in [3] where it was used 
in the context of relating the wave functions obtained in the Schrodinger quantization based 
on connection variables and the wave functions of the Wheeler-DeWitt quantization based 
on the metric variables. Since the two sets of basic variables differ significantly only when 
the spin connection is non-zero which happens when the intrinsic curvature is non-zero, non- 
trivial phases show up in this situation. The same transformation then appears whether one 
uses a loop quantization or a Schrodinger quantization in the connection formulation and 
compares with the Wheeler-DeWitt quantization. We also showed that the mechanism for 
the absence of singularities continues to hold, despite the phases, exactly as in the case 
of the Bianchi-I model [6]. Thus, as the main result of this paper we proved that for all 
diagonal Bianchi class A models there is a loop quantization which satisfies all the condi- 
tions for a good continuum limit as well as that of a singularity-free evolution. That such a 
quantization exists at all, and even one which provides a general form for all the models, is a 
non-trivial fact and gives further credence to the viability of loop methods and the physical 
applications obtained so far. 

In section HVl we carried out a preliminary exploration of the consequences of the quantized 
spin connection and the correspondingly modified potential. We focused on the Bianchi 
IX case and pointed out that unlike the classical potential which has infinite walls as the 
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singularity is approached, the quantum modified potential has finite walls. Since the infinite 
walls are responsible for the chaotic approach to the singularity, we expect that the finite 
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Is of the modified potential will significantly alter the approach to the classical singularity 
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Finally a remark about the relevance of the results is in order. Since the connection 
of the quantization discussed here to the full theory is weaker, it is legitimate to question 
what this quantization may have to do with the full theory. Firstly, we observe that the 
methods adopted as well as the two admissibility criteria work uniformly for all homogeneous 
models (including isotropic sub-models) which do not have a non-zero spin connection. The 
direct loss of connection to the full theory is only in the models with spin connection. The 
treatment of the spin connection enters in two distinct places, firstly in using the inverses 
of triad components and secondly in the modification of holonomies. The first one is quite 
natural and inherits the ambiguities in the definition of inverses The second one has 
additional ambiguities due to factor ordering. It is already clear that the first type of 
ambiguity is not fixed by the two admissibility criteria. Our emphasis in this work has been 
to demonstrate the existence of a uniform method of quantization which satisfies the two 
admissibility criteria and is singularity free. The uniqueness type of issues, while important, 
requires a systematic classification of possible ambiguities and is beyond the scope of the 
present work. Likewise the general question of how many of the results of a minisuperspace 
quantization can survive the full theory is also beyond the scope of this paper. 

While the spin connection is seen to require a careful treatment of the continuum ap- 
proximation and also to lead to a modified approach to the BKL singularity in the context 
of homogeneous models, its role in the full theory is more complicated and is an open issue. 
For instance, the phase factor in s which has the form, 

exp (z 7 - 1 ^ V IV) = exp (i 7 - V / d 3 xE?r a ) 

is not even well-defined in the full theory due to the transformation properties of the spin 
connection (which is no longer a covariant object). Since the spin connection can be made 
arbitrarily small locally, by choosing appropriate coordinates, one can expect that the dif- 
feomorphism constraint and its algebra with the Hamiltonian constraint will play a role in 
this issue. Requirements such as the local stability condition can also be expected to play a 
role in analyzing this issue, at least with a local version in the absence of a global internal 
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time. 
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